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A simplified physical model for calculating the onset temperature ratio and the frequency of a standing 
wave thermoacoustic engine (SWTE) in the time domain is built based on thermodynamic analysis. Coef¬ 
ficients of transient pressure drop and heat transfer are first deduced from linear thermoacoustic theory. 
By numerical computation, the evolutions of the pressure amplitude and the spectrum characteristics 
during the onset process are presented. Furthermore, the effects of stack spacing, charge pressure, and 
resonator length on the onset temperature ratio and the frequency are calculated. Relatively good agree¬ 
ment between the computational and the experimental results has been achieved, which validates the 
model for calculating the onset characteristics. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

A thermoacoustic engine converts heat into acoustic power by 
thermoacoustic effects. It has been attracting more and more 
attention for its high efficiency, simple configuration, high reliabil¬ 
ity, and capability to utilize low temperature thermal energy [1,2], 
Many efforts have been made to reveal the operating principle of 
thermoacoustic phenomena. Linear thermoacoustic theory, pro¬ 
posed by Rott [3], is now widely accepted. Only paying attention 
to the steady state and simplifying the analyzing process, Swift 
et al. converted the differential equations for thermoacoustics into 
algebraic equations by introducing complex notation [4,5], The 
time terms are eliminated in the linearized equations, which 
means that the equations cannot apply to the simulation in time 
domain directly. However, the time evolution process is indispens¬ 
able for understanding the mechanism of thermoacoustic onset 
process. 

In recent years, nonlinear thermoacoustic theories have been 
developed rapidly. The quasi-one-dimensional nonlinear models 
proposed by Yuan and Karpov [6,7] and the two-dimensional mod¬ 
el proposed by Hamilton [8] are capable of simulating the self- 
excited oscillation process of a SWTE in time domain. Meanwhile, 
computational fluid dynamics (CFD) has been used to simulate 
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thermoacoustic engines [9-12], It is worth noting that the 
transient characteristics of pressure drop and heat transfer are nec¬ 
essary for the above numerical simulations in time domain to 
achieve accurate results. However, there are few methods that sup¬ 
ply sufficient and reliable data for the oscillating flow in the stack 
or regenerator, due to the limited theoretical [13-16] and experi¬ 
mental [17-19] studies. Therefore, most nonlinear models have 
had to make use of the time-averaged data of pressure drop and 
heat transfer coefficients to calculate the energy and momentum 
exchange terms. 

A thermodynamic analysis method was proposed by de Waele 
[20] to reveal the onset conditions of a traveling wave thermoacou¬ 
stic engine, where all the components of the engine were con¬ 
verted analogously. To simplify the solution, an ideal regenerator 
was assumed, in which the energy equation was not required. 
The analytical expressions for the onset temperature and fre¬ 
quency were obtained. However, due to the intrinsically imperfect 
heat transfer in the stack of a SWTE, the method cannot be applied 
directly to SWTEs. 

In this study, all components of a SWTE are converted analo¬ 
gously to build a full model [21 ]. The energy equation is introduced 
to describe the imperfect heat transfer characteristics in the stack. 
Particularly, six coefficients were deduced from the linear ther¬ 
moacoustic theory to calculate the transient pressure drop and 
heat transfer in the stack. The whole process of building the model 
is presented in detail here. By solving the differential equation set 
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Nomenclature 



A 

area, m 2 

0q 

a coefficient for heat transfer, s/rad 

c 

sound velocity, m/s 

e T 

a coefficient for heat transfer, s/rad 

C P 

specific heat at constant pressure, J/(kg K) 

0v 

a coefficient for pressure drop, s/rad 

c 

flow conductance, m 3 /(Pa s) 

p 

density, kg/m 3 

D 

diameter, m; a coefficient for pressure drop, kg rad/ 

s K 

thermal penetration depth, m 


(m 4 s) 

s v 

viscous penetration depth, m 

f 

frequency, Hz 

T 

temperature ratio 

Sk 

spatial average thermal function 

CO 

angular frequency, rad/s 

fy 

spatial average viscous function 



G 

a coefficient for heat transfer, J/I< 

Prefix 


H 

a coefficient for heat transfer, J rad/(I< s) 

S 

deviation from average value 

L 

length, m 

A 

difference value 

rh 

mass flow rate, kg/s 

Subscript 

m 

mass, kg 

b 

buffer behind the resonator 

p 

pressure, Pa 

c 

cool end 

Pr 

Prandtl number 

g 

gas 

Q 

rate under which heat is transferred, W 

h 

hot end 

<3 

heat flow per unit area, W/m 2 

m 

mean 

R 

gas constant, J/( kg K) 

0 

Orifice 

T 

temperature, K 

R 

resonator 

u 

velocity, m/s 

s 

stack 

u 

volume flow rate, m 3 /s 

0 

time-average value for oscillation 

V 

volume, m 3 

l 

First order component for oscillation parameter 

y 

perpendicular distance from plate in the stack, m 


yo 

<l> 

y 

half of the stack spacing, m 
porosity of the stack 
specific heat ratio 

Special symbols 

() spatial average perpendicular to axial direction 

p 

dynamic viscosity, Pa s 




numerically, the onset processes of the SWTE are acquired. The on¬ 
set temperature ratio and oscillation frequency of the SWTE are ob¬ 
tained and further verified by experiments. 

2. Theoretical analysis of a SWTE 

Fig. 1 shows the schematic of a SWTE consisting of a hot tube, a 
hot heat exchanger (HHX), a stack, a cold heat exchanger (CHX), a 
resonator, and a buffer. Except for the resonator, the lengths of 
other components in the SWTE are usually much smaller than 
the acoustic wavelength. Thus, all the components except for the 
resonator are treated by lumped parameter method, in which com¬ 
pliance and resistance components are converted into volumes 
with uniform pressures and an orifice together with a buffer, 
respectively. Similar to de Waele’s treatment about the resonator 
[20], the resonator together with the buffer is converted into two 
volumes with uniform pressures and an in-between moving mass. 
Because of the intrinsic imperfection of the thermodynamic pro¬ 
cess occurring in the stack of the SWTE, coefficients deduced from 
linear thermoacoustic theory are used to calculate the transient 
pressure drop and the heat transfer in the stack. 

To simplify the solution and reflect the real conditions, theoret¬ 
ical analysis of the SWTE in this study is based on the following 
assumptions: 


(1) The pressure amplitude in the SWTE is much smaller than 
the charge pressure. 

(2) The solid temperature profile in the stack is linear and inde¬ 
pendent of time. 

(3) For the hot tube and the buffer, only the compliance effect is 
considered. For the resonator, only the inertance effect is 
considered. 

(4) The void volumes of the heat exchangers FIHX and CHX are 
ignored. 

(5) Except for the stack, the HHX, and the CHX, the other com¬ 
ponents are considered to be adiabatic. 

(6) There is no acoustic streaming occurring in the SWTE. 

(7) Working medium is assumed to be ideal gas. 

It should be pointed out that since the inertance effect of a short 
resonator is not dominant any more, omitting the compliance ef¬ 
fect and resistance of the resonator for a simplified solution may 
lead to a relatively large deviation when the resonator is short. 

As the resonator is so long that it cannot be considered as a 
compartment with uniform pressure, the components of the SWTE 
should be converted analogously to simplify the analysis. The con¬ 
verted SWTE is shown in Fig. 2. With the above assumption (3), the 
hot tube and the buffer can be regarded as compartments with uni¬ 
form pressure. The resonator is replaced by a solid piston with a 



Fig. 1 . Schematic of a SWTE. 
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mass m R . The orifice together with a buffer is a load of the SWTE, 
which may represent the dissipation in the resonator and the buf¬ 
fer behind the resonator. 

The lengths of the resonators at the two sides of the piston are 
L r and L' r , respectively. Let m R equals the mass of gas contained in 
the resonator in front of the piston, i.e. m R = p 0 L R A R . 


(a) If the resonator is connected to a buffer with infinite volume, 
the resonance angular frequency of the gas spring system 
combined with the piston and the compartments meets 
the following equation: 


co 


2 


yp 0 

PqLr 


( 1 ) 


For a tube with an open end, the resonance angular frequency 
m' is: 


c o' = 


TIC 

2 L a c 


( 2 ) 


where L ac is the length of the resonator. 

Making the resonance angular frequencies before and after the 
conversion the same, i.e. to = to’, then: 


U h 


1 dSPh 

Wh dt ’ 


( 8 ) 


where w h = ^ s» 22a. 

The parallel-plate type stack is used in this model for its high 
efficiency [5], According to the continuity equation of the stack: 


1 h-rh c = 


dm s 


> m< 


(\ dp s 
n-r— 


1 dT s 


U h - xU c = 


m s , o 

Ph 


'1 d5p s 
J>o dt 


1 dST s \ 


(9) 


( 10 ) 


The pressure drop of a thermoacoustic component derived by 
Swift [5] can be expressed as: 

dp, icop 0 icop 0 /A 

! X =-T^ u ' = ^r u '- (11) 


In time domain, the pressure drop along the stack can be ex¬ 
pressed as: 


SPh ~ s Ps = Dh 



U„, 


( 12 ) 


Lr = 


2 

71 


(3) 


(b) If the outlet of the resonator is closed, the resonance angular 
frequency of the gas spring system should meet the follow¬ 
ing equation: 


J3W1+1Y 

P 0 Lr \L r L r J 


(4) 


For a tube with close ends, the resonance angular frequency co' 
is: 


and Sp s - Sp c = D c 



U c , 


(13) 


where the coefficients D and 0 V indicate the pressure drop depen¬ 
dence of volume flow rate U and first-order derivation of U, 
respectively. 

Using the short stack approximation [5], it gives 


D|i(l + icod Vh ) 


1 ia >Pm h Ls 

2(1 -/vj V 


(14) 


271 2 L ac ' 

Making co = co' and combining with Eq. (3), then: 

L'r = \ Lr = 3 n Lac ' ( 6 ) 

The volume of the buffer after the conversion is: 

V„ = V' R + V’ b , (7) 

where V' b is the volume of the buffer before the conversion. 

Let T m = (T h + T c )/2, the pressure in the hot tube p h = p 0 + dp h , the 
pressure in the ambient end p c = Po + $Pa the pressure and the tem¬ 
perature in the stack p s = p 0 + Sp s and T s = T m + 3T S respectively. The 
pressure and the temperature in the stack are all space-averaged 
time-dependent parameters. 

As the characteristic dimension of the hot tube is usually much 
larger than the thermal penetration depth of the working gas, the 
thermodynamic processes in the hot tube can be taken as adia¬ 
batic. The volume flow rate U h is: 


andD c (l+ime Vc )=j^^_^ I (15) 

where A s = <j>A s is the gas area of the stack cross section. 

Separating the real and the imaginary parts, the coefficients D h , 
D c , 8 Vh , and 0 Vc can be calculated. 

From the equations derived by Swift [5], neglecting the small 
terms involving the conductivity and the heat capacity of the solid 
which will not affect much on the results, the temperature pertur¬ 
bation in the channel can be written as: 


7, = 


Jh 
Po c p . 
dr, 


l - 


cosh((l +i)y/d K ) 


cosh((l + i)2y 0 /5 K ) 

Pr cosh((l + i)y/S v ) 


1 dp | 
co 2 p 0 dx 


dx [Pr — 1 cosh((l +f)2y 0 /<5 v ) 
1 cosh((l +i)y/S K ) 


Pr - 1 cosh((l + i)2y 0 /5 K ) 


- 1 


(16) 


Using Eq. (16), the mean temperature over the channel width is: 



Fig. 2. Schematic of the SWTE after an analogy conversion. 
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(T,) = (1 -/*)-£r + 


1 


p 0 c p m 2 p 0 ( 1 - Pr) 
dT w dp^ 


X [Pr(l -/ v ) - (1 -f K [ 


(17) 


dx dx ' 

According to Eq. (16), the heat flux into the solid boundaries is: 


-qn = p 0 c p y 0 


i(of K 


Pi 


fv-f K 


dT n 


PoCp 0-/v)(l-Pr) dx 


u i 


(18) 


Combining Eqs. (17) and (18), the heat transferred from the 
boundaries is: 


0. = pc p V s 


imf K 


1 -u 


h 


1 dT w 
1 - Pr dx 


Pr - 1 


Pr 


l-/v 1 ~f K 


u 1 


From Eq. (19), the heat flow Qcan be expressed as: 

Q = Q(AT) + Q(u), 
where 

d 


Q(AT)=H(1+0 t -1(T s -T w ), 


and Q(u) =^G( 1 + 0, 


dt 


u. 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


Comparing with Eq. (19), we have: 

H(1 + icod T ) = pc p V s , 

1 Jk 

a„dC(, + ,W 0 )-aA(rr-1 +T i 7 -£ j . 

According to the energy equation of the stack: 

§(p s A s L s (c„T s + lu 2 )) +£[p s A s uL s (c p T s +±u 2 )] = 

§(c z ,m s T s + |m s u 2 ) +1 [( c p m s T s + \m s u 2 )L s ] = -[Q(AT) + Q(u)]. 

(25) 

Neglecting the small terms of the kinetic energy term and Q(u), Eq. 
(25) can be simplified as: 


d(c„m s T s ) 
dt 


= c p m h T h - c p m c T c - Q_(AT). 


Let m s T s = Eq. (26) can be expressed as: 

c P m h T h - c p m c T c + Q(AT). 

Dividing by ^ yields: 

= yRm h T h - yRm c T c - $-H( 1 + 0 T |)(r s - T w ) 

= y£p(u h -u c )-y=)H(\ + e T ^sT s . 


(26) 


(27) 


(28) 


dissipation. Thus, the thermal function in the resonator f KR is used 
to calculate the flow conductance: 


Co = C'Re(f KR ) 2 . (32) 

The motion equation of the piston in the resonator yields: 
d 2 x R 


mR - c[t 2 L = MdPc-SPb)- 


(33) 


The volume of the resonator compartment in front of the piston 
can be expressed as: V R = V R0 + A R x R . Eliminating x R yields: 


d 2 V R 

dt 2 


= a R (Sp c - Sp„), 


(34) 


where a R = -A. 

K m R 

Taking the thermodynamic process in the resonator as adia¬ 
batic, it gives: 


_dV R 1 dSp c 
R dt + w R dt 

where w R = ^~ y ^. 

According to the continuity equation: 

U c = U R + U„. 


(35) 


(36) 


Combining Eqs. (35) and (36) to eliminate U R yields: 




The dynamic equations of the buffer reveal: 

dV b = 1 dSp b 

dt w b dt ’ 


(37) 


(38) 


where w b = ^ ss For dV b = -dV R , Eq. (38) can be expressed as: 


dSp b 

dt 


w b 


dVR 

dt 


(39) 


For a SWTE with determined dimensions and operating param¬ 
eters, there are eight unknown variables in the equation set formed 
by Eqs. 8, 12, 13, 29, 30, 34, 37, and 39. Thus, the equation set is 
closed and can be solved numerically. 


3. Computational results 

Combining the above governing differential equations, a de¬ 
tailed case of a self-designed SWTE using nitrogen of 1.0 MPa as 
the working gas has been provided. The main dimensions of the 
SWTE are shown in Table 1 . 

To solve the governing equation set, let dV R /dt = F R . Substituting 
it into Eqs. 34, 37, and 39 yields: 


Combining Eqs. (10) and (28) yields: 

-gp = a^Uh + a 2 U c + a 3 ST s , (29) 

and ^ = b,U h + b 2 U c + b 3 5T s , (30) 

where a, = PoX^7Po + XilW, a 2 = -p 0 Zi(l’Po + XilTd a 3 = -(y 
- 1 bi =p 0 T m Xi(y - TJT h ), b 2 = PoTmX-i(T m /Tc - y), b 3 = 

Zi-l/[PoV s + (y-l)T m HO r ]. and & = (?-! )T 2 m 

H9 t /V s . 

The volume flow rate through the orifice U 0 is: 

U 0 = C 0 6p c . (31) 

The orifice, together with the buffer, represents the dissipations in 
the resonator and the buffer, including viscous and thermal 


dpR 

dt 


a R (Sp c -dp b ), 


(40) 


dSp^ 

dt 


= w b F R , 


( 41 ) 


and d ^=w R (U c -U 0 -F R ). (42) 

Eqs. (8), (12), (13), (29), (30), (40)-(42) form a complete set of 
eight first-order differential equations, and classical Runke-Kutta 
integration is used as the algorithmic method. An initial perturbance 
is needed to excite the oscillation at the beginning. Values of Sp c = 0, 
Sp h = 0, Sp s = 0, Sp b = 0, U h = 0, U c = 0, 5T S = 0, and F R = 1.0 x 10~ 5 m 3 /s 
are given as the initial values. The conductance coefficient C is set to 
be 1.0 x 10~ 4 m 3 /(Pa s). The cold end temperature ofCHX is set to be 
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Table 1 

Main dimensions of the self-designed SWTE. 


Components 

Dimensions 

Hot duct 

D/, = 59.5 mm, L b = 60 mm 

Stack 

D s = 59.5 mm, l s = 110 mm, y 0 = 0.3 mm, 0 = 0.4 

Resonator 

D r = 39 mm, L ac = 2600 mm 

Buffer 

D b = 85 mm, L b = 500 mm 


300 K. An initial angular frequency co is presumed to start the oscil¬ 
lation process and then the calculated oscillation angular frequency 
is used to get the final results. 

To calculate the onset temperature, the SWTE is simulated with 
different heating temperatures, increasing from 300 K step by step. 
When the pressure amplitude starts to grow, the heating temper¬ 
ature is taken as the onset temperature. Figs. 3 and 4 show the 
pressure evolution and the spectrum analysis in the hot tube with 
the temperature of the HHX of 450 K and 650 K, respectively. The 
given initial pressure amplitude damps gradually when T h is 
450 K, which indicates the heating temperature is lower than the 
required onset temperature. Whereas the pressure amplitude 
grows gradually when T h is 650 K, which means the heating tem¬ 
perature is higher than the required onset temperature. According 
to the spectrum analysis diagram, the oscillation frequencies for 
the two heating temperatures are 41.16 Hz and 41.32 Hz, 



(a) Pressure evolution diagram in hot tube 



(b) Spectrum diagram in hot tube 

Fig. 3. Pressure evolution and spectrum diagrams in hot tube for T h = 450 K and 
C = 1.0 x 10 4 . 



(a) Pressure evolution diagram in hot tube 



Fig. 4. Pressure evolution and spectrum diagrams in hot tube for T h = 650 K and 
C = 1.0 x 10 4 

respectively. The oscillation frequency of the converted resonance 
piston and the buffer is/ = ^/(wr + w b )a R /2n = 42.82 Hz, which is 
a little larger than the frequencies in simulations. This result indi¬ 
cates that the oscillation frequency of the SWTE system is mainly 
dependent on the resonant piston and buffer. The small difference 
between the oscillation frequencies of the SWTE system and the 
resonant piston and the buffer is due to the effects of the volumes 
of the hot tube and the stack. The operating frequency increases 
with the heating temperature because the physical properties of 
working gas in the hot tube and the stack change. As shown in 
the two spectrum diagrams, there is no obvious high order har¬ 
monic frequency, i.e., the energy is mostly focused in the funda¬ 
mental acoustic component. In fact, there is usually a second- 
order harmonic frequency observed experimentally [22,23], The 
difference is caused by replacing the resonator with a resonance 
piston in the model, which cannot show the harmonic 
characteristics. 

Onset temperature ratio is a critical parameter to determine the 
heating temperature needed for a thermoacoustic engine to start. 
Typically, lower onset temperature ratio indicates that it is more 
likely for a thermoacoustic engine to utilize lower grade thermal 
energy. For a SWTE operates depending on a nonideal heat transfer 
between the solid and working gas, the stack spacing is one of the 
most important parameters in the design of a standing wave ther¬ 
moacoustic engine. The effects of stack spacing on the onset tem¬ 
perature ratio of the SWTE are given in Fig. 5. As shown, the 


















DM. Sun et at./Applied Acoustics 81 (2014) 50-57 


55 


onset temperature ratio varies with the stack spacing for a given 
charge pressure. There exists an optimal y 0 at which the lowest on¬ 
set temperature ratio can be achieved corresponding to each 
charge pressure. The optimal y 0 decreases as the charge pressure 
increases under the present conditions. Furthermore, the onset 
temperature ratio varies with charge pressure for a given stack 
spacing. When the stack spacing is smaller than the optimal y 0 cor¬ 
responding to the highest charge pressure, i.e. 2.1 MPa in the pres¬ 
ent case, a higher charge pressure results in a lower onset 
temperature ratio, otherwise the opposite. Table 2 shows the opti¬ 
mal values of y 0 with different charge pressures. It is shown that 
the lowest onset temperature ratios for different charge pressure 
are the same (about 1.708) under the present conditions. The ratio 
yolS K for the lowest onset temperature ratio ranges from 1.28 to 
1.33, which indicates that the optimal stack spacing is about 2.6- 
2.7 times the thermal penetration depth. 


4. Experimental verification 

In order to verify the above model, a SWTE was constructed and 
tested. The onset temperatures and pressure oscillations were 
measured at various charge pressures for different resonator 
lengths. The main dimensions of the SWTE are listed in Table 1, ex¬ 
cept for the resonator length. Seven resonators of different lengths, 
i.e. 0.6 m, 0.8 m, 1.1 m, 1.6 m, 2.1 m, 2.6 m, and 3.1 m, are tested in 
the experiments. 

A calibrated piezoresistive pressure sensor with a total measur¬ 
ing error of about ±643 Pa was installed to measure pressure. The 
heating temperatures were measured by a calibrated NiCr-NiSi 
thermocouple with a maximum measuring error of ±3.5 K. Due 
to the difficulty in measuring the gas temperature directly, the so¬ 
lid temperature of the HHX is measured and used instead of the ac¬ 
tual gas temperature. The solid temperature is higher than the gas 
temperature. 

The effects of resonator length on the oscillation frequency with 
a charge pressure of 1.0 MPa are shown in Fig. 6. Both the calcu¬ 
lated and experimentally obtained frequencies decrease as the res¬ 
onator length increases. The calculated frequencies are larger than 
the experimental values because the compliance in the resonator is 
neglected. As stated before, the deviation between the calculated 
and experimental frequencies increases as the resonator length de¬ 
creases because of the relatively more remarkable compliance ef¬ 
fect for shorter resonators. 

Fig. 7 shows the calculated and experimental results of the ef¬ 
fects of charge pressure on the onset temperature ratio with 



Fig. 5. Effects of stack spacing y 0 on onset temperature ratio with different charge 
pressures (calculated results). 


Table 2 

Optimal parameters to reach the lowest onset temperature ratio. 


Po/MPa 

1.0 

1.2 

1.5 

1.8 

2.1 

Lowest onset temperature ratio 

1.708 

1.708 

1.708 

1.708 

1.708 

Optimal yo/mm 

0.225 

0.2 

0.18 

0.17 

0.15 

Optimal S K I mm 

0.172 

0.155 

0.139 

0.127 

0.117 

Optimal y 0 /<5 K 

1.31 

1.29 

1.29 

1.33 

1.28 



Length of resonator (m) 


Fig. 6. Effects of resonator length on frequency with the charge pressure of 1.0 MPa. 

different resonator length. The variation tendency of the calculated 
onset temperature ratio agrees well with what is obtained in 
experiments, especially when the length of the resonator is longer 
than 1.1 m. Both the calculated and the experimental results show 
that a lowest onset temperature can be achieved at a specific 
charge pressure for all the resonator lengths, which is defined as 
the optimal charge pressure here. The ranges of the optimal charge 
pressures of the simulations and experiments agree well. Similarly, 
there also exists an optimal resonator length to reach the lowest 
onset temperature for the SWTE, which is about 0.8 m in the calcu¬ 
lations but around 1.1 m in the experiments. 

The difference between the calculated and the experimental on¬ 
set temperature ratios are shown in Fig. 8. The corresponding 
charge pressures of the maximum deviation are given above the 
points. As shown, both the maximum and the average deviations 
are within 15.5% when the resonator is longer than 1.6 m. For a res¬ 
onator length of 3.1 m, the maximum deviation is only 7.4%, which 



Charge Pressure (MPa) 

Fig. 7. Effects of charge pressure on onset temperature ratio with different 
resonator length when the stack spacing y 0 is 0.3 mm. The experimental and the 
computational results are denoted by the filled and open symbols, respectively. 
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shows an accurate prediction of the model for long resonator 
SWTEs. The deviations increase as resonator length decreases, 
which is similar to the frequency deviations. For the shortest reso¬ 
nator, the maximum deviation is near 35%. The deviations between 
calculated and measured onset temperature ratios are mainly due 
to two reasons. Firstly, the measured onset temperature is the solid 
wall temperature of the HHX, which is higher than actual gas tem¬ 
perature. Secondly, neglecting the compliance effect of the resona¬ 
tor results in the large deviation when the resonator is relatively 
short. Furthermore, it is found in Fig. 8 that the maximum devia¬ 
tions for resonators longer than 1.6 m occur at relatively low 
charge pressures and resonators shorter than 1.6 m at relatively 
high charge pressures. Although the inertance effect of long reso¬ 
nators is dominant compared to the compliance effect, the superi¬ 
ority of long resonators in inertance effect decreases when charge 
pressure gets smaller, which causes larger deviations because only 
inertance effect of the resonator are considered in the model. On 
the other hand, as the length of the resonator decreases, the com¬ 
pliance effect becomes more and more remarkable. For example, 
when the length of the resonator is 1.1 m and the charge pressure 
is 1.0 MPa, the magnitudes of the compliance impedance and the 
inertance impedance of the resonator are 4.24 MPa s/m 3 and 
2.61 MPa s/m 3 , respectively. As the compliance impedance in¬ 
creases with the charge pressure, it is so remarkable for the short 
resonators at large charge pressures that neglecting it also causes 
larger deviations. 

In order to further verify the model, experiments done by other 
researchers are also included for comparison. Arnott et al. built a 
SWTE, whose structure and dimensions were reported in detail 
[24]. The cold end heat exchanger and resonator were cooled by 
water to maintain the temperature at 293.15 K. The onset fre¬ 
quency and the onset temperature ratio of the SWTE are calculated 
and compared with the experimental results presented in Ref. [24], 
as shown in Figs. 9 and 10. Fig. 9 shows that the frequency of the 
SWTE is almost independent of charge pressure. The calculated fre¬ 
quencies are larger than the experimental values by about 12%. In 
Arnott’s SWTE, the diameters of the hot tube and the resonator 
were almost the same, and the length of the resonator was about 
1.29 m. Omitting the effects of compliance of the resonator and 
inertance of the hot tube causes the deviation between the calcu¬ 
lated and experimental frequencies. Fig. 10 shows that the calcu¬ 
lated onset temperature ratio agrees well with that of 
experiments in tendency. There exists an optimal charge pressure 
to achieve the lowest onset temperature. The maximum difference 
between the calculated and experimental onset temperature ratio 
is about 11.9%. 



Resonator Length (m) 


Fig. 8. Effects of resonator length on deviation of onset temperature ratio. 



Charge Pressure (MPa) 


Fig. 9. Effects of charge pressure on onset frequency for Arnott’s SWTE. 



5. Conclusion 

This paper presents a simplified full physical model to simulate 
the onset characteristics of a SWTE based on thermodynamic anal¬ 
ysis. Coefficients deduced from linear thermoacoustic theory are 
used to calculate the transient pressure drop and heat transfer in 
the stack. The growth conditions of pressure oscillation, evolutions 
of pressure amplitude, and the spectrum characteristics during the 
onset process are obtained. The calculated frequency and the onset 
temperature ratio agree well with the experimental results espe¬ 
cially for long resonator SWTE. 
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